First-principles study of a;-Pu203 
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We systematically investigate the electronic structure, magnetic order, and valence states of a- 
PU2O3 (C-type) by using first- principles calculations. Q-PU2O3 can be constructed from Pu02 by 
removing 25% oxygen atoms. Our results show that the Pu 5/ orbitals are further localized after 
removing ordered oxygen atoms. This phenomenon is demonstrated by the combined fact that (i) 
the volume per unit cell expands 7% and (ii) the corresponding magnetic moments and valence states 
for Pu ions increase and decrease, respectively. According to the density of states and charge density 
distribution analysis, PUO2 is found to be more covalent than q:-Pu203, which is also because of the 
more localization of 5/ orbitals in the latter. The calculated lattice constants, bulk modulus, and 
electronic structures for PUO2 and a-Pu203 are consistent well with experimental observations. 

PACS numbers: 71.27.+a, 71.15.Mb, 71.20.-b 
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I. INTRODUCTION 



Plutonium-based materials have been extensively in- 
vestigated because of not only their great technological 
importance in the nuclear industry application but also 
their interesting rich physical properties from the basic 
theoretical viewpoint. Metallic Pu locates in the bound- 
ary of localized and delocalized 5/ electrons among the 
actinide metals and it has six different phase under differ- 
ent temperatures and pressures because of the complex 
character of the 5/ electrons [J, As for plutonium 
oxides, which are the only products when metallic plu- 
tonium is exposed in air, can store the surplus metallic 
plutonium The corrosion oxidation of metallic plu- 
tonium in air is very hot topic since it is a key problem 
for the protection and storage of plutonium-based nu- 
clear weapons. Furthermore, the thermodynamics and 
redox properties of plutonium oxides are also complex 
and interesting. When Pu is exposed to dry air at room 
temperature, the plutonium dioxide PUO2 layer is formed 
and a thin layer of plutonium oxide PU2O3 is followed on 
plutonium surfaces 0]. After several months or years, 
most of the Pu02 layer auto-reduces into PU2O3 layer. 
Especially, at 150-200°C, all the PUO2 layer auto-reduces 
into PU2O3 layer in minutes [1]. Usually, the mentioned 
plutonium sesquioxide is /3-PU2O3 in the hexagonal struc- 
ture (P3mi). However, another abnormal body cen- 
tered cubic form, plutonium sesquioxide PuOi.52, was 
also detected Ql- The ideal stoichiometric cubic pluto- 
nium sesquioxide is a-Pu203 (C-type PuOi.5), which are 
consisted of 32 Pu atoms and 48 O atoms per unit cell. 
a-Pu203 has not been prepared as a single-phase com- 
pound, since it is stable only below 300°C 5]. A mixture 
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FIG. 1: (left) unit cell of PUO2 with CaF2-type structure, 
(right) unit cell of a-Pu203 (body-centered cubic). The blue 
and red spheres denote Pu and O atoms, respectively. 



of cubic PuOi.52 and PuOi.gs can be obtained by partial 
reduction of PUO2 at high temperature and cooling to 
room temperature 

Naturally, the reduction of PUO2 into a-Pu203 is easy 
to happen because of the similarity of their cubic struc- 
tures. Pu02 crystalHzes in the CaF2 structure (see Fig. 
1) with the plutonium and oxygen atoms forming face- 
centered and simple cubic sublattices, respectively, and 
the cubic a-Pu203 can be obtained from Pu02 2x2x2 
supercell after removing ordered 25% oxygen atoms 0]. 
Actually, the structure of q;-Pu203 is body-centered cu- 
bic with space group /a3 (No. 206). Oxygen atoms occupy 
the 48 e sites, and plutonium atoms occupy the 24rf and 
8a sites [5] (see Fig. 1). Furthermore, if the 16c sites are 
also occupied by oxygen atoms, the resulting cell is the 
same as PUO2 2x2x2 supercell. Therefore, as mentioned 
above, the cubic a-Pu203 can be obtained by removing 
oxygen atoms located in the 16c (0.25,0.25,0.25) sites. 

In this work, we pay special attention to how the elec- 
tronic structure, magnetic order, and Pu valence states 
change during the Pu02-Pu203(a) reduction process. In 
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Pu02, Pu is in the ideal stable +4 oxidation state while 
in a-Pu203 the ideal oxidation state is +3. This change 
of valence states may illustrates that the 5/ electrons are 
more localized in a-Pu203. Similar to Ce02-Ce203 re- 
duction transition 6] , after creating one oxygen vacancy, 
two electrons left behind are condensed into localized /- 
level traps on two plutonium atoms, therefore, the va- 
lence for Pu changes from -t-4 to -1-3. Furthermore, due to 
the localization our results show that the volume expands 
7% during PUO2-PU2O3 isostructure reduction process. 

In order to describe more accurately the strong on- 
site Coulomb repulsion interaction among the plutonium 
5/ electrons, we use the generalized gradient approxima- 
tion {GGA)+U scheme. Our previous studies showed 
that the GGA-I- U approach can accurately describe the 
electronic structures and thermodynamic properties of 
Pu02 and /3-PU2O3 [8], which has motivated more the- 
oretical calculations in these two years [lol [Til . [T3 |. In 
the following, our calculated results demonstrate that 
the GGA-I- U correction can also successfully predict the 
ground state properties of cubic a-Pu203 and reliably de- 
scribe the Pu02-Pu203(a) reduction process. Our study 
shows that the system's volume expands during the re- 
duction process. Also, two well-resolved peaks in the den- 
sity of states (DOS) are observed and identified to origi- 
nate from the Pu 5/ and O 2p states in PUO2, which are 
consistent well with the recent photoemission measure- 
ments 0. In addition, our calculated bulk modulus for 
Pu02 is 180 GPa with Hubbard U^S.O eV, which agrees 
excellentl y w ith the recent refined experimental values of 
178 GPa jisj . However, in Ref. [ij], this value is largely 
overestimated to be 379 GPa. 



II. DETAILS OF CALCULATION 

Our first-principles calculations are based on the den- 
sity functional theory (DFT) and the Vienna ab initio 
simulation package (VASP) [15;] using the GGA for the 
exchange correlation potential [lB|. The electron and 
core interactions are included using the frozen-core pro- 
jected augmented wave (PAW) approach which combines 
the accuracy of augmented-plane-wave methods with the 
efficiency of the pseudo-potential approach [l3|. The Pu 
5/, 6s, 6p, 6d, and 7s as well as the oxygen 2s and 2p 
electrons are explicitly treated as valence electrons. The 
electron wave function is expanded in plane waves up to 
a cutoff energy of 500 eV. For the Brillouin zone integra- 
tion, the 2x2x2 Monkhorst-Pack sampling is adopted. 
The strong on-site Coulomb repulsion among the local- 
ized Pu 5/ electrons is described by using the formalism 
formulated by Dudarev et al. (3]. In this scheme, only 
the difference between the spherically averaged screened 
Coulomb energy U and the exchange energy J is impor- 
tant for the total LDA (GGA) energy functional. Thus, 
in the following we label them as one single effective pa- 
rameter U for brevity. In our calculation, we use J=0.75 
eV for the exchange energy, which is close to the values 



used in other previous work [1, [Tl| . 



III. RESULTS AND DISCUSSIONS 
A. Atomic and electronic structure of PUO2 

In order to investigate the electronic and structural 
properties of q;-Pu203 and how these properties change 
during the reduction process from Pu02 to a-Pu203, it 
is essential to first calculate the corresponding proper- 
ties of Pu02. Table I shows our calculated lattice con- 
stant Co and bulk modulus Bq and its pressure deriva- 
tive Bq within GGA-I- U scheme with different Hubbard 
U parameters. We performed the ferromagnetic (FM) 
and anti-ferromagnetic (AFM) coupling calculations to 
decide which is the ground state according to the total 
energies for each choice of the Hubbard U value. PUO2 
is known to be an insulator [l^ and some scattered ex- 
perimental data proved its ground state to be an AFM 
phase [1^. For bare GGA, i.e., U=Q eV, our calculated 
result predicts Pu02 to be a FM metal, similar to the 
conclusion of previous studies. This is contrary to ex- 
perimental observation. By increasing the amplitude of 
[/, as shown in Fig. 2, the total energy difference be- 
tween the two phases decreases; the AFM phase begins 
to be energetically preferred at U^l.h eV. We find that 
the convergence to the correct AFM phase is very rapid. 
Taking [7=3.0 eV, for example, the total energy (per 
formula unit) for the AFM phase is lower than for the 
FM phase by a large value of 0.4 eV. Note that in our 
calculation of Pu02 in this section, we have chosen two 
kinds of magnetic configurations. One is that the Pu 
magnetic moments fx are confined along the z axis in a 

simple H 1 — alternation of spins, while the other is the 

innerlayer alternation of spins. Our results show that the 
total energies for the two configurations are degenerate. 
In the ideal ionic limit, the formal charge for Pu ions in 
Pu02 is +4, corresponding to the formal population of 
for / orbitals. Our direct calculated magnetic moments 
/i for each plutonimn ion is 4.16 (in unit of Bohr magne- 
ton) mainly contributed by 4/ orbitals of 4.09, which is 
very close to the population of 

The equilibrium lattice constant ao and bulk modulus 
-Bo and its pressure derivative Bg are obtained by fitting 
the Murnaghan equation of state [l^]. At U =3.0 and 4.0 
eV, the present lattice constants ao are 5.458 and 5.477 
A for the AFM phase, respectively, which are in good 
agreement with the experimental value of 5.396 A[21|. 
As for the bulk modulus, our calculated results are 180 
and 184 GPa with U =3.0 and 4.0 eV, respectively, which 
agree well with the recent refined experimental value of 
179 GPa. However, in an early experiment this 
value is largely overestimated to be 379 GPa. Notice that 
the results of lattice constant and bulk modulus are 5.46 
A and 220 GPa, respectively, which are obtained using 
hybrid-density- functional calculations [2^ . The physi- 
cal insulating behavior of PUO2 with an experimentally 
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FIG. 2: The total energies for Pu02 in FM and AFM phases 
with different Hubbard U parameters. 



TABLE I: Calculated and experimental lattice constants 
ao(A) and bulk modules i?(GPa) for PU2O3 and PUO2 by 
the GGA+ U scheme. 
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4.0 


FM 


5.454 


184 
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AFM 


5.477 


184 


3.72 
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FM 


5.438 


183 
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AFM 
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3.95 







AFM 


5.398 


187 


3.67 
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FM 


10.91 
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10.92 


123 


4.25 




3.0 


FM 


11.14 


119 


4.05 




3.0 


AFM 


11.17 


122 


3.82 
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FM 


11.18 


120 


3.91 




4.0 


AFM 


11.20 


128 


3.31 




Expt. 
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comparable Mott gap of ~1.5 eV has as well been ob- 
tained with Hubbard J7 in a range of 3-4 eV. Therefore, 
our GGA-I- U results can provide a satisfactory descrip- 
tion of the atomic, mechanical, electronic, and magnetic 
(FM or AFM) structures of PUO2. This encourages us 
to investigate in the following the ground-state proper- 
ties of a-Pu203 by the same method, which turns to 
be remarkably effective as well. The present calculation 
of Pu02 is also necessary for theoretical access to the 
Pu02^Pu203(q;) redox energy, which will be discussed 
below. 



structure of q;-Pu203 is well confirmed by X-ray pow- 
der diffraction measurement [3, As mentioned above, 
after building the Pu02 2x2x2 supercell and removing 
25% oxygen atoms in special sites, q;-Pu203 can be ob- 
tained. Therefore, there are 32 plutonium atoms and 
48 oxygen atoms in the unit cell for q;-Pu203 (see Fig. 
1). The calculated structural and mechanical parame- 
ters within GGA-I- U scheme with different Hubbard U 
parameters for a-Pu203 are also collected in table I. We 
also performed the calculations in FM and AFM phases 
for a-Pu203. For the AFM calculations, we adopted two 
different spin arrangements of plutonium ions. One spin 
arrangement is to confine along the z axis in a simple 

intra- or interlayer H 1 — alternation of spins, as what 

we dealt with Pu02 in the last section. Another model 
is to use the same spin arrangement as in a-Mn2 03. 
Based on the neutron diffraction measurement, Regulski 
et al. [131 proposed a appropriate coUinear AFM order- 
ing, in which four basic magnetic cubic sublattices are 
constructed and the summation of magnetic moments of 
eight ions are assumed to be zero in each cube. We also 
adopt this AFM configuration as one additional candi- 
date to calculate the total energy of q;-Pu203 as a func- 
tion of U . The calculated lattice constant in the lat- 
ter AFM ordering are 10.92, 11.17, and 11.20 A for [/= 
0, 3.0, and 4.0 eV, and the corresponding bulk modu- 
lus are 123, 122, and 128 GPa, respectively. Note that 
the corresponding experimental lattice constant is in a 
range of 11.03 to 11.07 A[3, H. We find that similar to 
Pu02, at Hubbard U=Q eV, the result obtained by bare 
GGA predicts a-Pu203 to be a FM metal. Consider- 
ing that traditional DFT approach within the bare GGA 
scheme underestimates the strong on-site Coulomb re- 
pulsion of plutonium 5/ electrons and can not accurately 
describe the localization of 5/ electrons, therefore, we 
conclude this result obtained by bare GGA is not correct 
although no magnetic susceptibility or neutron powder 
diffraction data are available. Further calculations show 
that at Hubbard [/=3.0 eV and 4.0 eV, the ground state 
of a-Pu203 is an AFM insulator (see below). Notice 
that our direct calculations show that the q;-Pu203 in 
AFM phase is preferred to the latter spin arrangement as 
mentioned above. The corresponding energy differences 
Ai5=i?AFMi— -Bafm2 are 43 and 26 meV at C/=3.0 and 
4.0 eV per unit cell, respectively. As for the stability of 
AFM phase, the energy differences Ai^=i?AFM2— -^fm 
are —3.63 and —5.80 eV at U—i.Q and 4.0 eV per unit 
cell, respectively. 



Changes in the electronic properties during 
reduction process 



B. Atomic and electronic structure of a-Pu2 03 

a-Pu203 in body-centered cubic structure is the end 
product during the reduction action process. The crystal 



In Pu02 and q;-Pu203, the ideal oxidation states of 
Pu ions are -1-4 and -1-3, respectively. This valence differ- 
ence apparently shows that the population of 5/ changes 
and in depth demonstrates that the behavior of Pu 5/ 
is very complex. Directly, q;-Pu203 can be constructed 
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FIG. 3: The total DOS for AFM Pu02 and projected orbital- 
resolved partial DOS for Pu 6d, Pu 5/, and O 2p at Hubbard 
(a) U =0 eV and (b) U =3 eV, respectively. The Fermi level is 
set to zero. 
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FIG. 4: The total DOS for AFM a-Pu203 and projected 
orbital-resolved partial DOS for Pu 6d, Pu 5/, and O 2p at 
Hubbard (a) U =0 eV and (b) U —3 eV, respectively. The 
Fermi level is set to zero. 



from Pu02 by removing 25% oxygen atoms. After re- 
moving one oxygen atom, two additional electrons are 
left behind. These two electrons will occupy the low- 
est unoccupied states derived from Pu 5/ orbitals. This 
leads to the fact that the population of forms in a- 
PU2O3, while in Pu02. This clearly illustrates that 
Pu 5/ orbitals are more localized in a-Pu203 compared 



with Pu02. In the following, we systematically investi- 
gated how the electronic structures characterized by 5/ 
electrons change from Pu02 to a-Pu203. In order to 
make comparison, the electronic structure properties of 
Pu02 2x2x2 supercell containing 96 atoms are also cal- 
culated. We will show in detail how the volume of unit 
cell, density of states (DOS), charge distribution, valence 
state, and magnetic order change. All the following re- 
sults are obtained by GGA+U scheme with [/=3.0 eV 
except extra hints. 

(I) The calculated lattice constants for PUO2 and a- 
PU2O3 in their AFM phase are 5.458 and 11.17 A, leading 
to 7% volume expansion. This suggests that the 5/ elec- 
trons are more localized in Q!-Pu203. Due to the increas- 
ing localization, the interaction between Pu ions and the 
cohesion of the crystal decrease, therefore the lattice con- 
stant increases. The similar phenomenon has also been 
observed in the reduction process from Ce203 to Ce02 
with a 10% volume change 6] also because of the increas- 
ing localization of the / orbitals. 

(II) The total DOS for Pu02 (96 atoms) and a-Pu203 
(80 atoms) within GGA+ U scheme at different Hubbard 
U are shown in Figs. 3 and 4, respectively. In order 
to make a clear comparison, the orbital-resolved partial 
DOS for Pu 6d, Pu 5/, and O 2p orbitals are also plotted. 
As showed in Fig. 3(a), at U —0 eV, it is clearly that the 
bare GGA scheme predicts PUO2 to be a metal, which 
is contrary to the experimentally established insulating 
ground state [3. However, in Fig. 3(b), at Hubbard 
U=3.0 eV, one can see that the present result correctly 
predict PUO2 to be an AFM insulator. As for a-Pu203, 
at Hubbard U=0 eV showed in Fig. 4(a), the bare GGA 
approach predicts the FM metal ground state as we ex- 
pected, which is supposed to be not reasonable by con- 
sidering the necessity of the Hubbard U correction. At 
Hubbard U=3.0 eV, the calculated DOS showed in Fig. 
4(b) again suggests a-Pu203 to be in the insulating AFM 
ground state. Concerning the satisfactory description of 
the electronic structures for Pu02 and the fact that the 
5/ electrons in a-Pu203 are more localized than in Pu02, 
we predict a-Pu203 to be an insulator in the AFM state, 
although no experimental data can be obtained. 

As for the projected orbital-resolved partial DOS 
showed in Fig. 3(b) for Pu02, the occupied DOS is fea- 
tured by two well-resolved peaks. One near —1.8 eV is 
dominated by Pu 5/ character, while another one near 
—4.0 eV is mostly O 2p, which have observed in the re- 
cent photoemission measurement @. One can see that 
Pu 6d, Pu 5/, and O 2p states (especially between Pu 
5/ and O 2p) shows a significant hybridization covering 
from —5.7 to —0.4 eV. This hybridization is also respon- 
sible for the strong covalency of PUO2. However, for a- 
PU2O3, the Pu(5/)-0(2p) hybridization is much smaller 
due to the fact that Pu 5/ and O 2p occupied states are 
well separated as showed in Fig. 4(b). Furthermore, the 
covalency in q;-Pu203 is weaker than Pu02, which can 
also be demonstrated by the charge density distribution 
discussed below. We also note that in a-Pu203, the oc- 




(a) PuO, 



(b) Pu 



FIG. 5: Charge density in the (110) plane for (a) Pu02 and (b) Q-PU2O3, respectively. The contour lines are plotted from 
0.000 to 1.300 by the interval of 0.1083 e/k^. 



cupied 5/ peak around —1.3 eV is very narrow, indicating 
that the 5/ orbitals are very locaHzed, which is consistent 
with the conclusion in above part (I). 

(III) In order to obtain further understanding of the 
electronic structure and bonding properties, the contours 
of the charge densities in the (110) plane are also plotted 
for AFM PUO2 and AFM a-PusOa in Figs. 5(a) and 
5(b), respectively. Note that the contour lines are plotted 
from 0.000 to 1.300 by the interval of 0.1083 e/A^. For 
Pu02 in Fig. 5(a), it is evident that there are some 
closed contours existing between Pu and O atoms. This 
suggests that the covalent bonding character exists in 
Pu02 as discussed in above part (II). Compared with 
the charge density distribution for a-Pu203 showed in 
Fig. 5(b), one can see that the charge density at most 
bridges connecting Pu and O atoms are larger in Pu02 
than in ck-Pu2 03, this demonstrates that the covalency is 
much stronger in the former. Note that the distribution 
of charge density around Pu atoms is nearly spherical, 
while the distribution around O atoms is a little deformed 
towards their bonds. Therefore, it is easy to decide their 
ionic radii according the minimum value of charge density 
along the nearest Pu-0 bond. Then we can obtain the 
valence states for Pu ions discussed below. 

(IV) For PuOa and a-PuaOg, there are 12.463 and 
12.835 electrons around the plutonium ions with ionic 
radii of 1.287 and 1.319 A, therefore, the plutonium ions 
are presented as Pu^-54+ ^nd Pu^-^^+, which are close 
to the corresponding ideal valence states +4 and +3 for 
plutonium ions, respectively. 



(V) As for the magnetic order, we notice for Pu02 
(2x2x2 supercell) the second spin arrangement model 
is also more energetically favored with the energy dif- 
ference A£;=i;AFMi-£^AFM2 of 0.733 and 0.372 eV (per 
2x2x2 supercell) at Hubbard t/=3.0 and 4.0 eV, respec- 
tively. As for a-Pu203, the corresponding energy differ- 
ences A£^=i?AFMi--EAFM2 are 43 and 26 meV at U=3.0 
and 4.0 eV per unit cell, respectively. Our direct calcu- 
lations show that the magnetic moments are about 4.2 
and 5.0 fi^ per Pu ion in PUO2 and a-Pu203, which are 
mainly contributed by 5/ electrons of 4.1 and 5.0 fig and 
very close to the corresponding populations of and 
for plutonium ions, respectively. Notice that all the 
results of AFM phase listed in Figs. 3, 4, and 5 are ob- 
tained by the second spin arrangement model because it 
is energetically favored. 



D. Reduction reaction energy 

To further know the thermodynamic properties of the 
reduction process from Pu02 to q;-Pu203 via the reaction 

Pu02(2 x 2 X 2 supercell) - 8O2 = Pu203(a), 

we have calculated the reaction energy AE—Epu20a- 
£^Pu02+8£^02 at Hubbard U—3.0 eV. Note that Ep^^Oa, 
£'Pu02i and E02 are the total energies for PU2O3 (unit 
cell), Pu02 (2x2x2 supercell), and (spin-polarized) oxy- 
gen molecule, respectively. Considering the overestima- 
tion of the binding energy for O2 introduced by DFT, 
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we use the experimental cohesive energy value of 5.21 
eV for O2 [1^. Finally, the reaction energy is 55.04 eV 
(per a-Pu203 unit cell). This indicates that the reaction 
is endothermic, which is consistent with the fact that 
a-Pu203 can be experimentally obtained by partial re- 
duction of Pu02 at high temperature 

IV. SUMMARY 

In summary, we have studied the structural, electronic, 
and thermodynamics properties of PUO2 and a-Pu203 in 
their antiferromagnetic insulator states. We find that af- 
ter partial oxygen atoms removed from Pu02, the left 
electrons are localized into the Pu 5/ orbitals. As a con- 
sequence, several physical properties have changed cor- 
respondingly listed in following. (1) In PUO2, the ideal 
valence state for Pu is +A, while in a-Pu203, the valence 
state is -|-3. Our charge density integration around Pu ion 
sphere also confirms this conclusion. (2) The magnetic 



moments per Pu ion increase from 4 fig for Pu02 to 5 
fig for a-Pu203, which correspond to the populations of 
and respectively. (3) Due to the more localization 
of 5/ electrons, which do not contribute to the chemical 
bonding, thus the cohesion of a-Pu203 decreases. As a 
result, the volume expands about 7% according to our 
calculation. (4) Because of more localization of 5/ or- 
bitals, the occupied Pu 5/ and O 2p orbitals are well sepa- 
rated, the hybridization between them decreases. Finally 
the covalency of the Pu-0 bond in a-Pu203 is weakened. 
The calculated DOS and charge density distribution also 
demonstrate this. 
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